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Abstract 



An approach by Stephen is used to investigate the critical properties of 
randomly diluted resistor networks near the percolation threshold by means 
of renormalized field theory. We reformulate an existing field theory by Harris 
and Lubensky. By a decomposition of the principal Feynman diagrams we 
obtain a type of diagrams which again can be interpreted as resistor networks. 
This new interpretation provides for an alternative way of evaluating the 
Feynman diagrams for random resistor networks. We calculate the resistance 
crossover exponent up to second order in e = 6 — d, where d is the spatial 
dimension. Our result <j) = l+e/42+4e 2 /3087 verifies a previous calculation by 
Lubensky and Wang, which itself was based on the Potts-model formulation 
of the random resistor network. 

PACS numbers: 64.60.Ak, 64.60.Fr, 72.80.Ng, 05.70.Jk 
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I. INTRODUCTION 



Percolation is a leading paradigm for disorder. With regard to possible applications 
transport properties, e.g. electric conductance, are of particular interest in percolation the- 
ory. Electric transport on percolation clusters is commonly modelled by a random resistor 
network. In this model, bonds on a d- dimensional lattice are occupied with probability p 
and unoccupied with probability 1 — p. Each occupied bond has a finite nonzero conduc- 
tance u whereas unoccupied bonds have conductance zero. Suppose the system is near the 
percolation threshold, i.e. p is close to the critical concentration p c above which an infinite 
cluster exists. If one measures the resistance R(x, x') between two lattice sites x and x' 
known to be on the same cluster, one finds that the average over all possible configurations 
Mft obeys ~ |x — x'l*^, where v is the correlation length exponent defined by 

£ ~ {p ~~ Pc)~ v '■ The entire probability distribution for the resistance scales with the single 
exponent 0.^]. 

The theory of random resistor networks broke ground in the 70's [PHT2|j. Kasteleyn and 
Fortuin jl3| realized a connection between the random resistor network and the q — > limit 
of the g-state Potts-model. Stephen |TJ| found an approach connected to the xy-mode\. 
A Potts-model based formulation |]2],[T5|| enabled Lubensky and Wang fllfj to calculate the 
crossover exponent up to second order in e = 6 — d, d being the spatial dimension, as 
= 1 + e/42 + 4e 2 /3087. Nevertheless the approach by Stephen has been more fruitful |17 . 
For example, it has been employed to calculate for the random resistor network up to 



first order in e | l~7fl , several crossover exponents for a diluted network of Josephson junctions 
[18|, for a network with a singular distribution of resistances ||19fl , noise exponents for the 



random resistor network with fluctuating conductances |20[ as well as crossover and noise 
exponents for a randomly diluted network of nonlinear resistors f2~]| . However, we are not 
aware of any work up to now calculating based on this approach up to second order in e. 

The state of the art in the theory of random resistor networks dates back to the 80's. 
Since then progress has been reported only rarely, e.g. by Fourcade and Tremblay [E2 . The 
complexity that the field theory of random resistor networks has reached might be a reason 
for this. Here we reformulate a field theory by Harris and Lubensky JT7J based on Stephen in 
a way, that we believe is less complex and more intuitive. We present in detail a calculation 
of up to second order in e. We hope this fosters further calculations of this type which 
might appear in the future. 



II. THE MODEL 

In this section we line out the essentials of the derivation of a field theory for random 
resistor networks due to Stephen |14| and Harris and Lubensky [17]. We will provide the 



reader with indispensable background and clarify certain points. 



A. Kirchhoff 's equations 



Consider a ci-dimensional lattice. Each bond is occupied by a resistor of conductance 
cr with probability p or unoccupied with probability 1 — p. Moreover each lattice site is 
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connected to ground by a capacitor. For convenience all capacitors are taken to have the 
value unity. Kirchhoff 's law applied to site i reads 

Qi = E (ft - Vi) +Ii = xuoVi , (2.1) 
j 

where Qi and Vi are the charge and the potential at site % and u is the corresponding 
frequency, Ij is an externally imposed current and the sum runs over all nearest neighbors. 
In order to obtain a solution for the voltages Eq. ( |2.1|) may be cast into matrix form, 

gY_ = I, (2.2) 

by setting S^j = —o%j and S iti = iu + J2j For w ^ the inverse of g is well defined 
and we can write 

V = g^I ■ (2.3) 

The limit u —>■ requires some cautiousness. To characterize clusters C we define vectors 
e(C) with ej(C) = if % ^ C and e$(C) = NiC)' 1 ^ 2 if i G C, where N(C) denotes the number 
of sites belonging to C. The vectors e are eigenvectors of S_: 

g e(C) = iue{C) . (2.4) 

Thus the inverse of S_ may be written as 

Si.} — + > (2-5) 

with S_ 1 summarizing the part of S _1 that stays finite in the limit uo — > 0, ^ = A^(C) -1 
if i,j G C and = if z and j are not located on the same cluster. From Eq. ( |2.5|) we 
conclude that singularities due to zero modes can be regularized as far as the infinite cluster 
is concerned by taking the thermodynamic limit before uo — > 0. In the remainder of this 
article we focus on vanishing u unless stated otherwise. 

Suppose a current I is put into a cluster at site x and taken out at site x' . Those sites 
connected to both x and x' by two mutually non intersecting paths are constituting the 
backbone between x and x' . The current at a site i belonging to the backbone may be 
written as 



and the voltage at site i is simply 

The difference in voltage between the source and the sink is 



V x - V x , 



S~ l + S~ l , — 2S~\ 



(2.6) 



(2.7) 



(2.8) 



which states that the resistance R(x, x') between sites x and x' reads 



3 



R(x, x ) — S X)X + S x } x i 2>S XjX > ■ (2-9) 
Different cases may be distinguished: If x and x' lay on the same cluster then 

R(x^ x ) = S x x ~\~ S ' / x i x' ' ^2. 10) 

If x and x' are not connected then S~L = and hence 



R(x, x ) — — (z x>x + z x i ;X i) + (Sj.^ + >x , . 



(2.11) 



We see that R(x, x') — > oo for c<j — > if x and are not on the same cluster. In the limit 
of widely separated sites one has S~l, — > 0. Thus the resistance Roo between two infinitely 
separated points on the infinite cluster is 



R r 



(2.12) 



with the last equality holding for homogenous systems. 

Now we turn to the power P = I (V x — V x ') dissipated on the backbone. It may be 
written as 



P = R(x, x')~ l (V x - Vx>f = E °k,i (Vk - Vf = P ({V}) 



(2.13) 



k,l 



with the sum running over all nearest neighbor pairs on the backbone and {V^} denoting 
the corresponding set of voltages. In terms of P Kirchhoff's equation (|2.1| ) arises as a 
consequence of the variation principle 



_d_ 



-P({V}) + Y.hV k 







(2.14) 



Equivalent to Eq. ( 2.13 ) the power can be expressed in terms of the currents as 

P = R(x,x')P = J2^Uli- (2-15) 



k.l 



Suppose the backbone contains closed loops as sub-networks with currents {-f^} circulating 
independently around these loops. Then the current flowing through a certain bond is not 
only a function of / but also of the set of loop currents: 



(2.16) 



Conservation of charge holds for every ramification of the backbone and this gives rise to 
another variation principle: 







dP l ) 



p({/ w },/)=o. 



(2.17) 



Eq. (|2.17|) may be used to eliminate the loop currents and thus provides us with a method 
to determine the total resistance of the backbone via Eq. ( |2.15|) . 
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B. Replica formalism 



From the discussion above it is evident that our task is in principle to invert a random 
matrix. This inversion can be generated by gaussian integration. However, we are interested 
in the average resistance (R(x, x'))c an d hence the average over all possible realizations of the 
diluted configuration for fixed p remains to be performed. This can be achieved by employing 
replica technique [23]. The network is replicated D-fold: V_ — > V_ — {v}- l \ ■ ■ ■ ,Y^ D ^j. The 
replication procedure leeds to the extended generating function 

|z-°|nnc^(-^+a-£)| =(ex P (-ix-r 1 A)) c; , (2.18) 

with the power 

P = V ■ £ V = J2 Vl a) S itj v} a) (2. 19) 



and where X-V_ = Y,i, a K u ■ The normalization constant Z is adjusted in the usual way: 

Z= J]JdV i exp(-~V-iV) . (2.20) 

The choice Aj = A (5i, x — 8i tX >) provides us with a generating function for the average resis- 
tance: 

'^(x)^_ x (x')) rep = /exp I --R(x,x') ) ) , (2.21) 




where 



V> x (z) = exp (iX ■ V x ) ,A^0 (2.22) 




and (...)re P denotes the average over all replicated voltage and lattice configurations. An 
expression similar to Eq. ([2.21 ) holds for the resistance to infinity: 

(2.23) 

As we shall see later on, it is useful to define an indicator function that is unity if x and x' 
lay on the same cluster and zero otherwise. Consider the limit o — > oo. From the discussion 
in subsection [II A| it is clear that R (x, x') = if x and x' are connected and R (x, x') = oo 
else. This suggests to define the indicator function in terms of the average {ip^(x)ip_^(x')) 
over all replicated voltage configurations for a given realization of the diluted lattice as 

X (x, x') = \im(iP x (x)iP_ x (x')) . (2.24) 

Similarly 
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X (x) = fim (<i[> x (x)) (2.25) 

indicates if x is located on the infinite cluster. 

Some remarks should be made at this point. Since infinite voltage drops between different 
clusters may occur it is not guaranteed that Z stays finite, i.e. the limit lim^^o Z D is not 
well defined. This problem can be regularized by switching to voltage variables 9 taking 
discrete values on a D-dimensional torus. The voltages are discretized by setting 9 = A9k, 
where A9 = 9m/M is the gap between successive voltages, 9 M is a voltage cutoff, A; is a D- 
dimensional and M a positive integer. The components of k are restricted to — M < kS a > < 
M and periodic boundary conditions are generated by equating k^ a ' = k^mod(2M). The 
continuum may be restored by taking 9m — ► oo and A9 — > 0. By setting 9m = 9oM, M = m 2 
and respectively A9 = 9 /m the two limits can be taken simultaneously via m — > oo. Note 
that the limit D — > has to be taken before any other limit. Since the voltages and A are 
conjugated variables, A is affected by the discretization as well: 

A = AXt , AXA9 = n , (2.26) 

where I is a D-dimensional integer taking the same values as k. This choice guarantees that 
the completeness and orthogonality relations 

£ exp (i\-9)= 5 x g mod(2MAA) (2.27a) 

and 

exp (i\-6) = % ) o mod(2M Ae) ( 2 - 27b ) 

do hold. Eq. (|2.27|) provides us with the Fourier transform 

* f (x) = (2M)- D exp (iA ■ 9) tP x (x) = 8 §A - {2M)~ D (2.28) 

with the condition 

E¥ j; ) = - ( 2 - 29 ) 



Note that $ is nothing else than a Potts-spin EM with q = (2M) D states 



C. Field theoretic Hamiltonian 

We proceed with the evaluation of Eq. (|2.18| ). Carrying out the average over the diluted 
lattice configurations provides us with the weight exp(— H rcp ) of the average (...) rep : 

H iep = - In /exp (~P) \ = ^ % ■ % ~ £ ln ( ex P ("^m (S " *i) ' (6 " ^ )) ) ■ ( 2 - 30 ) 
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By dropping a constant term A^ln(l — p) with Nb being the number of bonds in the 
undiluted lattice we obtain 

#rep = Y, k (3) + E K ft ~ > ( 2 - 31 ) 

where 



and 




(2.32) 



(2.33) 



The Fourier transform of X, 




(2.34) 



can be Taylor expanded as 

oo 

^ (A) = r + 5> p (A 2 )" , (2.35) 

J3=l 

with r and u> p ~ a~ p being expansion coefficients. Since K decays exponentially in a9 2 
the series ( |2.35| ) may be terminated after the quadratic term for large a. Further on we 
omit factors (2M)~ D which become unity in the limit D — > 0. We define the discrete nabla 
operator Vg* through 

-E V AW • V^(x) = E AVx(s#_x(z) (2-36) 

and obtain 

K (A,-) = r - W A e - , (2.37) 

where w — W\. Similarly the right hand side of Eq. ( p.32| ) may be viewed as the leading 
term of an expansion 

oo 

h(e) = Y.KifY (2-38) 

with h\ = iuj/2. 

Now we can set up a field theoretic Hamiltonian TC in compliance with the symmetries 
of the model. The average over the configurations of the diluted lattice renders the model 
symmetric under spatial translations and rotations. Another feature of the model is that 
it is local since only nearest neighbors enter in Eq. (|2.1| ). In the limit of perfect transport 
(a — ► oo) and in the absence of external fields (h p = 0) the model is invariant against 



7 



permutations of all q = (2M) states of the Potts-spins $g(x). If one allows w p 7^ this 
S q symmetry is lost. The model remains gauge invariant under a shift of the voltages by an 
arbitrary finite potential as can be inferred from Eq. ( |2.1| ). This symmetry corresponds to 
translational invariance in replica space: &ff{x) «-> &q_q q (%)■ Moreover we deduce from the 
quadratic form of the power in Eq. (|2.19| ) that the model possesses of an 0(D) symmetry 



in replica space. Additionally admitting h p 7^ results in breaking of the gauge invariance. 

We proceed with the usual coarse graining step and replace the Potts-spins $ff(x) by 
the order parameter field if fx, (Pj . By constructing all possible invariants of the symmetries 

discussed above from J2gf ( x , &f an d gradients thereof the following Hamiltonian in spirit 
of the Landau-Ginzburg- Wilson functional is obtained: 

U = |^x^|^(x,^^(A^)^(x,^ + ^(x^y + /i(x,^^(x,^)} , (2.39) 

6 

where terms of higher order in the fields have been neglected since they turn out to be 
irrelevant. The coarse grained kernel and external field must resemble the features of the 
original K and h. This gives 

« = /^z (*> *y + f ( v ^ ( x ' + \ ( v ^ ( x ' ®)) 2 

e 

+ ^(x,^) 3 + ^(x^-)} , (2.40) 

where r, w and uo are now coarse grained analogues of the original coefficients. Note that 7i 
reduces to the usual (2M) D -states Potts-model Hamiltonian by setting w = h = as one 
one retrieves purely geometrical percolation in the limit of vanishing h and w (er — > 00) . 

It is worth pointing out that the problem of calculating the moments of the resistance 
distribution has been reshuffled: = ^x(x, x')i? /c (x, x ')) c ma y be obtained by taking the 
fc-th derivative of 



A 2 

lim (^ x (x)^__ x (x)\ = (x(x, x')) c - — (x(x, x')i?(x, x')) c + . . . 




2 

(2.41) 



with respect to A 2 at A 2 = 0. 



D. Scaling properties 

We conclude this section with a scaling analyses of the Hamiltonian (|2.39|) . Let P denote 
the set of parameters {w p , h p } and b some scaling factor for the voltage variable: 9 b9. 
By substitution of (p (x, (P) = <p' (x, b6j into the Hamiltonian we get 



S 



n 



<pf (x, be) , p] = [d d xJ2\ (x, be) k (a, Ag) iff (x, be) + |y (x, w) 



+ h(e)(p' (x, 60) 



^2.42) 



Renaming the scaled voltage variables 6 1 ' = be leads to 



' (x, e') ,p]=Jd d xj:\ \p' (x, e') k (a, 6 2 A e .) tf/ (x, 0') + |y (x, e') 3 

+ h (b- l ff) ip' (x, 6') | . (2.43) 

Clearly a scaling of the voltage variable results in a scaling of the voltage cutoff: 9m — > 60 a/ • 
However, by taking the limit -D — > and then m — > oo the dependance of the theory on the 
cutoff drops out. No ultraviolet divergencies are encountered in integrations over voltage 
variables and hence the voltages are no origin of anomalous scaling. We can identify 0' and 
e and thus 



n 



if (x, be) , p] = n [p (x, 9) , p' 



(2.44) 



where P' = {b 2p w p ,b 2p h p }. We can conclude the following implication of Eq. ( [2.44 ) on 
correlation functions 

G N ({x, 9} ; t, {w p , h p }) = Jvtptp(x x ,fy...<p (x N , 9 N ) exp (-H [p (x,8),P\) , (2.45) 
where Dip indicates an integration over the set of variables j<£> (x, 9) j for all x and 9: 

G N ({x, 9} ■ r, {w p , h p }) = G N ({x, bO} ; r, {b 2p w p , b~ 2p h p }) . (2.46) 
The two point correlation function G2 is the Fourier transform of (■0^(x) , 0_^(x)) aj . We 



deduce from Eq. ( (2.41| ) that 

\ 2 M R ((x, x') ; r, {w p , h p }) = (b~ 1 \) 2 M R ((x, x') ; r, {b 2p w p , b~ 2p h p ) 



H 



(2.47) 



The freedom to choose 6 has not been exploited yet. With the choice b 2 = w 1 the previous 
scaling relation turns into 



M R ((x, x') ; r, {w p , h p }) = wM R ((x, x') ; r, 0, w p 6, 



f2.48) 



For p > 1 coupling constants W7 P only appear as w p /w p . The associated exponent p0 — P 
is of order p — 1 as mean field analysis shows. This indicates that terms of order w p X 2p for 
p > 1 give rise to corrections to scaling. We keep only the leading terms w\ 2 and hi = iu/2. 
By virtue of w ~ a -1 the resulting scaling relation reads 



m h((x,x');t,ct,^) =a l f ((x, x') ; r, w/«t) 
where / is a scaling function of cu/cr. 



(2.49) 
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III. TWO LOOP CALCULATION AND RENORMALIZATION 



A. Diagrammatic expansion 



To perform the renormalization program we start out with a dimensional analysis. A 
trivial consequence of the fact that the Hamiltonian ( j2.4U| ) must be dimensionless is that 
the involved quantities have naive dimensions x ~ wX 2 ~ /i 2 , ip ~ /i^ -2 ^ 2 , r ~ /i 2 and 
g ~ fi ( - 6 ~ d ^ 2 , where /i is a convenient inverse length scale. The positive dimension of the 
coupling constant g shows that it is relevant for d < d c = 6. 

The principal elements contributing to the diagrammatic expansion are easily gathered 
from the Hamiltonian, namely the vertex —g and the propagator 

1 ~ ^A,0 = 1 ^3 - 

p 2 + r + wX 2 p 2 + t + wX 2 P 2 + r ' 

which is displayed in Fig. |l|. 

The superficial degree of divergence 5 in any one-particle-irreducible diagram composed 
of these elements is 5 = dL — 2P, where L denotes the number of loops and P the number 
of propagators. The topologic relations 3V = E + 2P and L = P — V + 1, with E being the 
number of external legs and V being the number of vertices, lead for d = d c to 5 = 6 — 2E. 
Therefore the only superficially divergent vertex functions are T 2 and T 3 (see Fig. |^). 



B. Feynman diagrams as resistor networks 



We learn from Eq. (|3.1| ) that the principal propagator (bold) decomposes into two prop- 
agators. One of them is carrying replica currents and we refer to it as conducting. The other 
one is not carrying A's and we call it insulating. This decomposition of the bold propagator 
allows for a schematic decomposition of bold diagrams into sums of diagrams consisting of 
conducting and insulating propagators. 

There is an important feature of the decomposition scheme that we want to point out at 
the instance of the diagram displayed in Fig. |3|. The diagram reads 

g 4 r 1 11 1 1 

~2 Y ^ q r + k 2 + w\ 2 r + k 2 r + k 2 r + (k + q) 2 + wk 2 t + q 2 + wk 2 ' ^ ' ' 

where J kq is an abbreviation for (27r) _d / d d kd d q. Schwinger parametrization leads to 

(3.2) = ^ E / r II ds i ex P {~ T E s i - (*i + s 2 + s 3 ) k 2 - s A (k + q) 2 - s 5 q 2 ) 

x exp (— (s 4 + s 5 ) wk 2 — siwX 2 ^j . (3.3) 
The sum over k factorizes 

E ex P (- («4 + s 5 ) wk 2 ) = I E ex P (~ ( s 4 + s 5 ) wk 2 ^ j (3.4) 
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and hence becomes unity in the limit D — > 0. What we encounter here in our example is a 
particular instance of the diagrammatic rule formulated in Fig. ^. 

We apply the decomposition scheme to all one and two loop diagrams. The result is 
displayed in Figs. |5|-f| 

From the decomposition a new interpretation of the Feynman diagrams emerges. They 
may be viewed as resistor networks themselves with conducting propagators corresponding 
to conductors and insulating propagators to open bonds. The Schwinger parameters Si 
correspond to resistances a" 1 and the replica variables — zAj to currents. As conservation 
of charge holds for every ramification in a resistor network, the A's are conserved in each 
vertex. The A-dependent part of a diagram can be expressed in terms of its power P: 



exp (-w s *^il = ex P ( wP (% {^})) 



(3.5) 



with Aj = Aj (X, {ft}), where (appart from a factor — i) A is an external current and {k} 
denotes the set of independent loop currents. 

The new interpretation suggests an alternative way of computing the Feynman diagrams. 
To evaluate sums over independent loop currents 

£exp( W P (%{£})) (3.6) 

{k} 

one can employ the saddle point method that is exact in our case since the power is quadratic 
in the currents. Note that the saddle point equation is nothing else than the variation 
principle stated in Eq. (p,17|) . Thus solving the saddle point equations is equivalent to 



determining the total resistance R ({s,}) of a diagram and the saddle point evaluation of 
(pi) yields 



exp (-i? ({s,}) wA 2 ) . (3.7) 

After a completion of squares in the momenta the momentum integrations are straightfor- 
ward. Thereafter all diagrams are of the form 

/ (p 2 , A 2 ) = I P (p 2 ) - I W (p 2 ) w\ 2 + . . . 
where D (p 2 , {sj}) is a usual integrand of the 3 -theory. 



C. Renormalization 



We use dimensional regularization ]25| to compute the various diagrams obtained by 
decomposition. Appendix [A] lines out these calculations in terms of examples. As the result 
of an e-expansion up to second order in e _1 we obtain for the superficially divergent parts 
of the vertex functions 
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T 2 (p, A; r, g, w) = r + p 2 + w\* - g^r^ { (l + r + ip 2 + \wX 2 ) 

d G 2 r/9 45 \ /ll 7 \ „ /65 169 , r , . 

+ £ 4 ^f - + — e r + — + e p 2 + — + e wX 2 \ 3.9a 

e 2 IV4 16/ V36 432 / H V36 432 ; ' 



and 



T 3 ({0, X};r,g,o)=g- 2g^r^ 2 + g^r^ (H + , ( 3 . 9 b) 

where G e = (47r) _d/,2 r(l + e/2) with T denoting the Gamma-function. We included the 
convergent term —g 2 r 1 ~ e / 2 G e in Eq. ( |3.9a| ) since it is important in the calculation of the two 
loop contribution to the corresponding Z-factor. The vertex functions ( p.9|) are apart from 
terms proportional w\ 2 identical to those of the Potts-model. 

The e-poles are compensated by minimal subtraction. We employ the following renor- 
malization scheme: 

ip ^ ^ = Z 1/2 tp , r^r = Z- x Z t t , (3.10a) 
w w = Z - l Z w w , g - g = Z^/ 2 Z^G^u 1 / 2 ^/ 2 . (3.10b) 

In minimal subtraction the Z-factors have to be determined such that they solely cancel 
the e-poles. We find 



z = 


1 • 


lu 
67 


37 u 2 11 u 2 _ N 
- 432 e + 3 6e2 + ° (M3) ' 


(3.11a) 


Z T = 


1 + 


u 
e 


47m 2 9m 2 ... „ 


(3.11b) 


Z w = 


1 + 


5 u 
67 


319 m 2 65 u 2 , nl ^ 
" 432 , + 36 £ » +0(u) ' 


(3.11c) 


Z u 


1 + 


a 
4- 

e 


59 M 2 n 2 „, . „. 

-127 + 11 ? + ° (M ) ' 


(3.11d) 



with Z, Z T and Z n being the known Potts-model Z-factors. 



IV. RENORMALIZATION GROUP EQUATION AND SCALING 

The unrenormalized theory has to be independent of the length scale introduced by 
renormalization, i.e. 

d 

^g-G N ({p,X};r,g,w) =0 (4.1) 

for all N. Eq. (|4.1|) translates via the Wilson-functions 

du d In r 

C(«) = ^, 7 («) = ^1bZ (4.2b) 
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(the bare quantities are kept fixed while taking the derivatives) into the Gell-Mann-Low 
renormalization group equation 



d n d d d N ' 

Vtt- + Ptt + tk— + w(— + —7 
o\i ou or aw 2 



G N ({x, X} ;t,u,w,h) = . 



(4.3) 



The particular form of the Wilson-functions can be extracted from the renormalization 
scheme and the Z-factors. We find 



P(u) = u(-u 



7 671 



72 



u 2 -e 



k (u) = -u u + O (u 3 ) , 

v ' 6 108 V ; ' 

2 47 

C(u) = u --u 2 + 0(u*) , 

3 3d 

1 37 

1 (U) = __ M+ _ M 2 + ( M 3) _ 



(4.4a) 
(4.4b) 
(4.4c) 
(4.4d) 



The renormalization group equation is solved by the method of characteristics. The charac- 
teristics read 







with 


m 


= V , 


(4.5a) 


du 


= (3{u{l)) 


with 


u(l) 


= u , 


(4.5b) 


l m lnr 


= «(«(«)) 


with 


f(l) 


= T , 


(4.5c) 


l—lnw- 
al 


= C(«(0) 


with 


w(l) 


= W , 


(4.5d) 




= 7(«(0) 


with 


Z(l) 


= 1 . 


(4.5e) 



Solving the first one is trivial. For the remaining characteristics fixed point solutions are 
determined. The fixed point condition (3{u*) = leads to the infrared stable fixed point 



* 2 671 



2 1 Oie 3 ) . 
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We obtain as fixed point solution of the renormalization group equation 
G N ({x, X} ; r, u, w, fi) = F* N/2 G N ({x, X} ; rl K \ u*, wl c , fd) , 



(4.6) 



(4.7) 



where 7* = 7 (u*), n* = k (u*) and C* = C ( u *)- 

To get a scaling relation for the vertex functions a dimensional analysis remains to be 
performed. It yields 

(4.8) 



G N ({x, X} ; r, u, w, //) = n (d ' 2)N/2 G N ({/xx, iT l w l,2 \) ; fi' 2 r, u,l,l) . 
From Eq. ( [4.7| ) and Eq. ( [4.8| ) we derive the scaling relation 
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G N ({x, X} ; r, u, w, //) = ^ d ~ 2+ ^ N / 2 G N ({Zx, / c/2_1 A} ; r^r, u*, w, //) (4.9) 
with the well known critical exponents for percolation 

77 = 7* = -— e-— e 2 + 0(e 3 ) (4.10) 

' r 21 9261 V ; K 1 

and 

„ = (2 - k*)- 1 = - + -e + -^-e 2 + O (e3) . (4.11) 

v ; 2 84 37044 v ; v ; 

(p remains to be determined. A glance at Eq. ( |2.41j ) shows that 

M£ ~ A" 2 . (4.12) 
The scaling properties of A can be deduced from Eq. (|4.9|) and hence 

M\ ~ |x-x'| 2_c * . (4.13) 



Thus we finally obtain 



u(2-C) = 1 + — e + — ^e 2 + (e 3 ) (4.14) 
v s ; 42 3087 V ; V ; 



in conformity with the result by Harris and Lubensky. 



V. RATIONAL APPROXIMATION 

Since the exact value of <fi is known to be unity not only in d > 6 dimensions but also in 
one dimension it suggests itself to do a rational approximation. The feature <f>{d = 1) = 1 is 
incorporated by rewriting (p as 

<P = 1 + e (— + —e - ^—e 2 + O (e 3 )) (5.1) 
v V42 3087 154350 V V V ; 

being identical up to second order to the e-expansion result. In Fig. |9] we compare the 
analytic result to numerical estimates for d = 2 by Grassberger [32] and d = 3 by Gingold 



and Lobb ||33|| . The extrapolated value <f>(d = 2) = 1.04 deviates from the simulation result 
0.9825 ± 0.008 by roughly 5%. In d — 3 one obtains <fi(d = 3) = 1.05 compared to the 
numerical estimate 1.117 ± 0.019. Here the deviation is about 6%. 

The good agreement should be taken with caution. Due to the rich structure of the 
exponent rj in the percolation problem, the exponent ip = ^ = <p / v — 2 + rj might be better 
suited for comparison to simulations. Rational approximation of ip yields 

/ 5 187 334 2 , r n , 

t/> = e eH e 2 + C(e 3 ) . (5.2) 

V V 21 3087 15435 y 1 y 1 



This is compared to the simulations in Fig. 10. The agreement for d = 3 is reasonable. As 



expected the discrepancy increases for d decreasing. For d = 2 the analytic result looks 
unrealistic. The structure of ip appears to be too rich to be reproduced at low dimensions 
by a series of a few terms. 
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VI. CONCLUSIONS AND OUTLOOK 



We have presented a study of randomly diluted resistor networks based on an approach 
by Stephen. The motivation has been twofold. Firstly we wanted to verify a result for the 
resistance crossover exponent </> obtained by Lubensky and Wang using a different approach. 
Our result for is in absolute agreement to that by Lubensky and Wang. Secondly we wanted 
to simplify the theory of random resistor networks. We demonstrated how a decomposition 
scheme leeds to a new interpretation of the involved Feynman diagrams: they may themselves 
be viewed as resistor networks. The new interpretation of the Feynman diagrams greatly 
improves the handling of the calculations. 

We are positive that our formulation will foster future investigations in transport on 
percolation clusters. As we will report shortly in a separate publication, it proved to be 
advantageous for nonlinear random resistor networks for which V ~ I r . In the limit r — ► +0 
the resistance between two points becomes essentially equal to the length of the shortest 
paths between these points. Our result for the exponent for this so-called chemical distance 
d min = 2 - e/6 - [937/588 + 45/49 (In 2 - 9/10 In 3)] (e/6) 2 + O (e*) verifies a previous calcu- 
lation by Janssen |26j . The limit r — > oo is related to the red (singly connected) bonds. Our 
two loop calculation gives unity for the corresponding exponent in accordance with results by 
Blumenfeld and Aharony ]27J and de Arcangelis et al. [f^|] . Moreover our formulation enabled 
us to calculate the percolation backbone exponent Db to third order in e by considering the 
limit r -> -1. We find D B = 2 + e/21 - 172e 2 /9261 + 2 (-74639 + 22680C (3)) e 3 /4084101, 
where ( denotes the Riemann zeta function. 
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APPENDIX A: EVALUATION OF DIAGRAMS 

In this Appendix we want to sketch the computation of the Feynman diagrams. The 
evaluation of contributions to T3 is straightforward. By virtue of 5 = one can set external 
momenta and wX's equal to zero. Hence all diagrams resulting from one bold diagram are 
giving the same contribution and the decomposition reduces to a mere factor. An example 
is given in Fig. [TI|. The diagrams obtained in this fashion are the usual ones found in the 
literature on the Potts-model p9f and can be evaluated by standard procedures |3(| . 



Now we turn to IV Since T2 is superficially divergent with degree S = 2 a Taylor 
expansion up to first order in p 2 and wX 2 is sufficient. Contributions of order zero and pro- 
portional to p 2 are again standard. Thus we restrict ourself to demonstrate the computation 
of contributions to T 2 proportional to wX 2 . 

We revisit our example of section |T| and start with Eq. ( |3.3|) . As we have concluded the 
sum over k merely gives a factor unity. After a completion of squares in the momenta the 
momentum integrations are straightforward. We get 
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(3.3) 



U ds 



exp 



(-tE*=i 



2 (47r) d Jo i=1 [( S4 + S5 ) ( Sl + S2 + s 3 + S4 ) - s | 



21^/2 



exp (siwX" 



At this stage it is useful to perform a change of variables: S4 — > ii, s 5 — ► t 2 , Si 
s 2 - ^ *32/ and s 3 — > £ 3 (1 — x — y). In these variables the integral reads 



(Al) 

• t 3 x, 



(Al) 



2 (4tt) 



dtidt 2 dt 3 / cfx 



exp(-r (ti+ta + t 3 )) 
Mi + hh + tit 2 



d/2 



exp (-t 3 xwX 2 ) . (A2) 



After expansion for small wX 2 the integrations with respect to x and y are easily carried 
out. We omit the term of order zero and obtain 



(A2) 



g A wX 2 
" 2~(4^) d Jo 



dt X dt 2 dt 3 l eX P(-^^+^ + ^)) t 3 



6 [t 3 ti + t 3 t 2 + ht 2 



d/2 



3 ' 



(A3) 



Eq. 



) can be expressed in terms of a parameter integral 

1 



M 1 (a, 6, c) 



p,q (a + p 2 ) (b + q 2 ) (c + (p + q) 



(47T) 



dt\dt 2 dt- 



exp (- (ati + 5t 2 + ct 3 )) 
[t 3 ti + £3^2 + ^1^2 



d/2 



(A4) 



by taking partial derivatives with respect to a, b or c at a = 6 = c = r. This parameter 
integral was introduced by Breuer and Janssen ||31|| . They find in dimensional regularization 



M 1 (a, 6, c) 



6e 



fl 25\ / 



fl 3-e + b 3~e + ^ 



3-e 



21\ r 



7 + 4/ 



a 2 " e (b + c) + (a + c) + c^ e (a + 6) 



3abc 



(A5) 



In terms of this parameter integral Eq. 

(A3) = -wX 2 ^ 
and is easily evaluated yielding 

(A6) = -w\ 



reads 



1 d^M 1 



6 <9c 3 



(A6) 



a=b=c=T - 



„4 r~<2 
2 9 ^ e 

2 e 



-r ^< - 



1 3 
6~e ~ 8 



(A7) 



As a second example we take the rightmost diagram in the first line of Fig. ||. This 
diagram corresponds to the resistor network in Fig. 12. The total resistance of this network 
is 



Hence the saddle point evaluation of 



(■Sl + S 2 ) (S3 + S 4 ) 

(si + s 2 + s 3 + s 4 ) 



(A8) 
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2 -V 



kq J 



X 



JJ rfsj exp -r ^ Si - (si + s 3 ) k 2 - s 5 q 2 

i=l V i=l / 

exp ^— (s 2 + s 4 ) (k + q) 2 - (si + s 2 ) wk 2 - (s 3 + s 4 ) (/? + X) ^ 



(A9) 



gives 



(A9) 



g 4 wX 2 [°° ^ ^ ( Sl + s 2 ) (s 3 + s 4 



2 (47r) d 7o 



i= i \ s i + S2 + s 3 + Si) 
exp (-r Ef=i Si 



x 



(A10) 



[s 5 (si + s 2 + s 3 + s 4 ) + (s 2 + s 4 ) (si + s 3 )] d/2 

where we have already carried out the momentum integrations and the expansion for small 
wX 2 . The change of variables s 2 — > tiX, s 4 — > £i(l — x), si — > t 2 y, s 3 — > £ 2 (1 — y) and S5 — > £ 3 
recasts the integral into 

g 4 wX 2 roc 



(A10) 



dt\dt 2 dt% I dxdy t\t 2 ^ ' 



x 



2 (47r) d 7o Jo 
(t l x + t 2 y) (t l (l-x)+t 2 (l-y)) 



(All) 



Carrying out the integrations with respect to x and y yields 

Since a similar structure emerges in several diagrams we introduce an additional parameter 
integral 



M (a,b,c) 



1 



dtidt 2 dU 



tit 2 exp (- (ati + bt 2 + e£ 3 )) 



(A13) 



(4vr) d y — — tl + t2 [ t 3 tl + f3t2 + tlt2 ]"/2 
We calculate M 2 (a, 6, c) in dimensional regularization which yields 

M 2 (a, 6, c) = S I (I + ( a *- + ^ + * Q + ^ + a 6 + 1 (a + 6) c} . (A14) 
Now we can evaluate ( |A12| ) by taking derivatives of M 2 (a, 6, c): 



(A12) 
We finally obtain 



-wX 



,# 4 f 1 <9 2 M 



2 3 da 2 



+ 



1 d 2 M 



a=b=c=r 



2 <9a<% 



a=b=c=T - 



(A15) 



?2 <7 4 G 2 _jl 1 
wA - — -t e { h - 



2 e 



3e 6 



(A15) 



(A16) 



The remaining diagrams two loop diagrams contributing to T 2 can be treated in a similar 
fashion. However the two diagrams in Fig. containing only light propagators are calculated 
most conveniently by observing that momenta and replica variables play exactly the same 
role in these diagrams (see Fig. [L3|). Appendix [FJ gives an overview of the two loop diagrams 
contributing to T 2 in terms of M 1 and M 2 . 
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APPENDIX B: THE DIAGRAMS IN TERMS OF PARAMETER INTEGRALS 



Here we list our results for the diagrams contributing to T 2 , as far as not given in 
Appendix For convenience we use the notation 



M W = (i _ i) !(j _ !)!(/ - i)! da^d^dd-^ 1,2 ( °' 6 ' C) 



(Bl) 



a=b=c=T 



and 



p (r + p 



:b2) 



Those parts of the diagrams proportional to w\ 2 are displayed in Fig. TA. The remaining 
parts can be inferred from literature on the Potts-model. 
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FIGURES 



FIG. | 

The principal propagator (bold) decomposes into two propagators. One of them (light) is 
carrying currents and we refer to it as conducting. The other one (dashed) is not carrying 
A's and we call it insulating. 

FIG. | 

The superficially divergent vertex functions T 2 and T 3 . Terms of 0(3 loop) have been ne- 
glected. 

FIG. | 

A diagram encountered by decomposing the contributions to T 2 . 
FIG. | 

An important feature of the decomposition of the bold diagrams. By applying the de- 
composition scheme one may obtain sub-diagrams which are connected to the rest of their 
diagram solely by insulating legs. The hatched blob on the left hand side stands for such a 
sub-diagram. For each closed loop of conducting propagators a sum over an independent k 
has to be performed. In the limit D — > these sums merely produce factors unity. Thus 
all conducting propagators can be replaced by insulating ones and there are no summations 
necessary in the sub-diagram. 

FIG. | 

Decomposition of the one loop diagram contributing to T 2 - 
FIG. | 

Decomposition of the two loop diagrams contributing to T 2 . 
FIG. 

Decomposition of the one loop diagram contributing to T 3 . 
FIG. | 

Decomposition of the two loop diagrams contributing to T 3 . 
FIG. | 

Dependence of the exponent <ft on dimensionality. The e-expansion up first (diamonds) and 
second order (squares) as well as the rational approximation (triangles) are compared to 
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numerical results (circles) by Grassberger and Gingold and Lobb. They determined the 
exponent t for the conductivity S, X ~ \p — p c |', by simulations, t is related to via 
= t — (d — 2)u. In d = 2 u is known exactly f34j|3"5[| : v = 4/3. For d = 3 we use a Monte 
Carlo result by Ziff and Stell @: v = 0.875 ± 0.008. 



FIG. |10| 

Dependence of the exponent tp — (t — 7)/// + 2 — <i on dimensionality. As in Fig. |9] diamonds 
and squares refer to the e-expansion, triangles to the rational approximation and circles to 
numerical results by Grassberger and Gingold and Lobb. For the exponent 7 governing the 
mass of finite clusters we use 7 = 43/18 |j34[ , [3~5| for d = 2 and 7 = 1.795 ± 0.005 [0] for 
d = 3. 



FIG. |TT 

Since T3 is superficially divergent with 5 = it is sufficient to evaluate the diagrams at 
zero external momenta and currents. Thus all diagrams resulting from one bold three leg 
diagram are giving the same contribution. The decomposition results in a mere factor, being 
identical to the tensor contraction in the Potts-model. 



FIG. 11 



The resistor network corresponding to the rightmost diagram in the first line of Fig. || The 
conducting propagators are interpreted as conductors whereas the insulating propagators 
are interpreted as open bonds. The Schwinger parameters Si correspond to resistances. 



FIG. IT31 



A convenient way to extract the contributions proportional to w\ 2 from diagrams not con- 
taining insulating propagators. 



FIG. Ill 



Listing of contributions to the diagrammatic expansion proportional to w\ 2 . The right 
hand sides remain to be multiplied by a factor —wX 2 G 2 /e which we dropped for notational 
simplicity. 



22 




k,A 

FIG. 3. 




FIG. 5. 
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FIG. 8. 
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FIG. 9. 




FIG. 10. 




FIG. 11. 




FIG. 12. 




p=0 
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FIG. 14. 



